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Abstract. We use multipolar expansions to find the force on a gold coated 
sphere above a gold substrate; we study both an empty gold shell and a gold 
coated polystyrene sphere. We find four characteristic separation ranges. In the 
first region, which for the empty gold shell occurs for distances, d, smaller than 
the thickness of the coating, the result agrees with that on a solid gold sphere and 
varies as d~ 2 ; for larger separations there is a region where the force behaves as 
if the coating is strictly two dimensional and varies as d~ 5 / 2 ; in the third region 
the dependence is more unspecific; in the forth region when d is larger than the 
radius, the force varies as d~ 4 . For homogeneous objects of more general shapes 
we introduce a numerical method based on the solution of an integral equation 
for the electric field over a system of objects with arbitrary shapes. We study 
the effect of shape and orientation on the van der Waals interaction between an 
object and a substrate and between two objects. 



1. Introduction 

In his seminar paper 1] Casimir calculated the energy between two parallel 
perfect conductor half spaces as the change of zero-point energy of the classical 
electromagnetic field, 



where (z) are the frequencies of the electromagnetic field modes of the plates 
separated by a distance z. This procedure has been employed by van Kampen [5] to 
study vdW (van der Waals) forces of dielectric plates in the non-retarded case. Later 
Gerlach [3] interpreted the retarded Casimir forces in terms of interacting surface 
plasmons of the plates. Most recently the problem of Casimir energy between two 
dielectric plates has been revisited and it has been shown that it can be obtained as the 
sum of contributions from the interacting surface plasmons and from the propagating 
modes of a cavity formed by the parallel plates [HOIS]. Apart from the original results 
for parallel plates and slabs, a limited number of results exists for other geometries and 
combination of materials. The most well known result is the so called proximity force 
approximation (PFA) [7]. In that model the interaction between objects of various 
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shapes is calculated in terms of the energy between planar interfaces. In the case of a 
spherical object of radius R, above a planar substrate, the force F becomes 
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where z is the closest distance and E p (h) is the interaction energy per unit area 
between two planar interfaces a distance h apart. The first factor of the result is what 
one usually means with the PFA. The second factor is a correction which we neglect in 
this work. In the case of vdW forces other analytical results exist for highly symmetric 
geometries. One of them is the following general expression for the interaction energy 
between two cylinders, spheres or half spaces at small separations 
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where n = for spheres, n — 1 for cylinders, n = 2 for half spaces and A12 is the so 
called Hamakcr constant given by 
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It has been shown that for larger separations, contrary to this formula, the energy is 
not symmetric with the interchange of materials between the objects [8]. The results 
of (|3|) in fact coincide with those obtained with the usual PFA approach for those 
geometries, in the nonretarded case which we are concerned with here. 

Another set of results use additivity of pair interactions between the atoms of two 
objets to obtain the behavior of the force at close and at long distances. In spite of 
the well known fact that the Casimir force is not an additive interaction, there is a 
general consensus that these results give the correct distance dependence of the forces. 
One example is that for two cylinders; for parallel cylinders the vdW energy goes like 
~ d~ 15 when d — > 0, but when they are crossed, a dependence ~ d~ 2 is expected, 
which is the same as between two spheres. 

A spectral representation formalism has been employed [9] to obtain the resonance 
frequencies of a sphere above a substrate from which one obtains the vdW energy. In 
the spectral representation formalism, the optical response of a two face system is 
expressed as the sum of resonant terms from which it is possible to explicitly obtain 
its resonance frequencies [10]. Those resonances correspond to the surface plasmons 
of the system. Furthermore, in that formalism the effects of geometry and material 
are clearly separated; this makes it possible to deduce the importance of the geometry 
in the resonances of the system, and consequently in its vdW energy. 

In this work we use a multipolar method similar to that in [5] to calculate the 
vdW interaction between a substrate at arbitrary separation from a spherical object 
with an internal structure in the form of a sphere with a coat of a different material. 
We must mention that there are no results for vdW interactions between finite objects 
with a structure formed by parts of different materials. We also introduce a formalism 
for the determination of the vdW energy of a system of objects of arbitrary shape, 
separation and orientation. We present results for the vdW forces between cubes, and 
finite cylinders above substrates, and develop calculations of lateral and rotational 
forces between objects. 
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2. Formalism 



2.1. Multipolar method 

First we consider the system shown in the inset of figure [TJ a). A system of four 
materials bounded by two concentric spheres of radii R — 8 and R with the center of 
the spheres at a distance R + z from the surface of a substrate. The material of the 
interior of the sphere with radius R — 6 has a dielectric function £j(w); the coat or 
shell bounded by R — S and R contains a material with dielectric function e c (cu); the 
substrate has a dielectric function e s (oj). The full sphere is immersed in an ambient 
with dielectric function e a (u>). It has been shown [91 111) that the multipolar moment 
Qi m of the surface charge induced by the fluctuations of the electromagnetic field 
with multipolar component V^ c of the stochastic vacuum field, satisfies the system 
of linear equations 



+ fcA?Az) 



h m = -V?Z, |m| </,/', (5) 

V 

where a; is the polarizability of the sphere (independent of m), f c = 
[e s (u) — c a (ui)] / [e s (ui) + e a (u))], and A™, are the terms, of the matrix that relate 
the multipolar functions between the different systems of spherical coordinates of the 
sphere and its image, that vanish when z — ► [TT1 [S] . For a homogeneous sphere ai 
is given by 
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where a is the radius of the sphere and u(iv) = [1 — e sp (uj) /e a (w)] _1 is the so called 
spectral variable. The function s sp (uj) is the dielectric function of the sphere material. 
From ^ and ^ one can obtain a spectral representation for Q; m , 
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are the so called spectral functions, where U™, 



is the orthogonal matrix that satisfies 



i,i' 



\z)8 s 



(8) 



Here, H™, (z) — n®5u> + J ' C y/c~iy/ci> 'A™, (z) is a real symmetric matrix. The resonance 



frequencies are given by 



and the vdW energy is obtained from ©. 



Note that the matrix H™, (and consequently n s ) depends only on the geometrical 
properties of the system and the dielectric properties of the substrate. In the case of 
a coated sphere the polarizability can be obtained using the boundary conditions of 
the electric potential as 

c (uj)} - [n, - u a M] (1 - 5/R) 2l+1 
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[nj - u ic (uj)} [m - u ca {uj)] +nj(l - nj)(l - 5/R) 2l+1 ' 
with u xy (uj) = 1/ [1 — e x (u>) / s v (lj)]. Now it is neither possible to find a spectral 
representation of the system nor to obtain its resonances explicitly. One must rely 
on the argument principle to calculate the energy of the system, as follows. The 
determinant of ([5]) becomes zero at the resonance frequencies of the system, 

' -6 lM + fMz)] =0. (10) 
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It is possible to apply the argument theorem to the multipolar problem [3] and using 
the fact that when z — > oo the resonances are those of the isolated coated sphere to 
find that ([T]) is given by 
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This expression can also be used in the case of experimental dielectric functions or 
complex valued model dielectric functions for the sphere and substrate. 



2.2. Integral equation method 

The multipolar method above gives very high precision results over a wide range of 
distances. However to deal with particles of other shapes the multipolar method 
becomes much more complex and expensive numerically. On the other hand for finite 
objects of arbitrary geometry it is not possible to define one unique way of calculating 
the energy in the PFA approach. For spheres and cylinders one considers small planes 
that interact with the opposite in the direction parallel to the line that joins the centers 
of the objects. However for more complex geometries it is not clear how to define the 
interacting planes. R Fuchs has developed a spectral representation formalism for the 
calculation of the optical response, in the non-retarded limit, of homogeneous objects 
of arbitrary shapes [12] . We use his basic formalism for the formulation of a model to 
obtain the vdW interaction energy for systems of objects of arbitrary geometry. 

Fuchs has shown that the surface charge density, a, satisfies the integral equation 
defined over the surface S of the object [T2"] . 

E~(r, W ) • n(r) = (r) - / ^'^ a (r) , (12) 

e(w)-1 J s |r — r'| 

where E 1,ac is the stochastic vacuum field that excites the particle and n is the 
surface normal pointing outwards from the surface. The complete deduction of (|12p 
is independent of the form and connectivity of the surface [12]. We use this integral 
equation when £ is the combination of the surfaces of two finite arbitrary objects. The 
simplest numerical solution of (fT2jl is obtained, as in |12j . by dividing the surface £ 
in small fragments of area Asi and surface normal and with homogeneous surface 
charge a i: centered at the points r,. As a result, one obtains the linear system of 
equations, 

27T— -Oi - > Rijo-i = fi, Ri,j = (1 - dij) ■ ^ , (13) 

and fi = E°(ri) • n(r^). A spectral representation can be obtained from solving IT3l) . 

X ( l4 ) 

^ u(u) - n s 

where n s = (1 — m s /27r)/2, m s are the proper values of the matrix Rij, and C] ,J are 
obtained from the matrix that diagonalizes Rij [T^]. The resonance frequencies are 
given by u(uj s ) = n s and like in the case of the multipolar spectral representation of 
a sphere above a substrate, the vdW energy could be calculated with ([5]). 
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Figure 1. (a) The radius times the force on a gold shell of radius R and thickness 
8 at distance d above a gold substrate as function of d/R. See the text for details, 
(b) The interaction energy for a gold coated polystyrene sphere of radius R and 
coat thickness S at distance d above a gold substrate as function of d/R. The coat 
thickness takes on the values 0.0, 0.001, 0.01, 0.05, 0.1, 0.3, 0.6 and 1, in units 
of R, counted from below. The thick solid curve is the PFA result based on two 
gold half spaces. 

3. Results and Discussion 

3.1. Free standing spherical shells and gold coated polystyrene spheres 

For the calculations in this part we first construct the matrix of the linear system, ([5]) , 
for a given z / R. For the dielectric function of gold and polystyrene we use those in 
[T3"] and [2] , respectively. We calculate numerically the determinant and perform the 
integration of An increasing maximum number L max (I, I' < L max ) of multipoles 
is considered until the convergency of the energy is reached. 

In figure [lja) we show the results for a free standing gold shell above a solid 
gold substrate. As long as one may neglect retardation effects the figure is universal, 
independent of the size of the sphere. The thick solid curve is the PFA result when 
E p (z), of (O, is the interaction energy between two half spaces; here the force varies as 
dr 2 . We see from the figure that the full result from using the multipolar method, thin 
solid curves, merge with the thick solid curve for small enough separations. This only 
happens inside the figure for the thickest of the coatings. For the other two examples 
it happens for smaller separations. For larger separations there is a region where the 
result follows the PFA result when E p (z) is the interaction energy between one half 
space and one strictly two dimensional film, dotted curves; here the force varies as 
d~ 5 / 2 , i.e., has a fractional power dependence. The dotted curves are obtained as 

[T32CL6] 

Rx F = 2ttR 2 E p (d) re Q.lhhQ\Jnh 2 e 2 /me^fo T jR~l {d/Rf /2 . (15) 

The circles are the PFA result when E p (z) is the interaction energy between a half 
space and a film of thickness S. Note that the last two calculations (dotted lines and 
circles) are extensions of the PFA model; in the original model one used E p {z) of two 
half spaces. The dashed curves are the result from only including dipolar interactions. 

Measurements of the vdW and Casimir interactions often involve spheres since 
in this case the difficult problem of alignment is obsolete. In most cases these spheres 
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Figure 2. vdW forces between gold prisms the distance d above a gold substrate. 
The thin solid curve is the result for a gold cube of side L. The dashed (dash- 
dotted) curve is the result for a standing up (lying down) prism of length 2L. The 
thick solid curve is the PFA result for the cube, or standing up prism, based on 
two gold half spaces. 



are not solid but coated; one assumes that the coated spheres may represent solid 
spheres. As we have seen, this assumption is only valid for relatively small distances. 
In figure D^b) we show the interaction energy for a gold coated polystyrene sphere 
above a gold substrate. This way of displaying the results also provides universal 
curves, independent of the radius of the sphere. The radius, R, is the value for the full 
sphere including the coating. We see that only the spheres with very thick coatings 
behave as solid gold spheres. 

3.2. Interaction between objects of arbitrary geometry 

3.2.1. Objects above substrates: For the calculation of the vdW energy of these 
systems we define a discretization of the surface of the object in small areas specified 
by their positions, normal vectors and sizes. With that information we construct the 
matrix Ri j in <| 1 3j) and obtain the eigenvalues m s numerically and from them the 
resonance frequencies. Finally, the energy is calculated with (O. Using a maximum of 
1000 discrete elements per object we reached down to distances of 2% of its maximum 
diameter, within a numerical error of less than 2%, as estimated from comparisons 
with calculations of [HI E] ■ 

In figure[2]we show calculations for different prisms above substrates, of the same 
material. The thin solid line shows the result for a cube of side L at a distance d; the 
dashed line shows the result for a prism of height 2L and square base with side length 
L; the thick solid line shows the result of ([3]). We note that when the cube and the 
square prism are very close to the substrate the force equals the base area times the 
force per unit area between two half spaces. At long distance the force is dominated by 
dipolar interactions which are proportional to the volume. The dash-dotted line shows 
the result when the prism is lying down, i.e., one of the rectangular sides is parallel to 
the substrate. So for small separations the force on the cube and standing up prism 
are equal while it is twice the size on the lying down prism; for large separations the 
forces on the standing up and lying down prisms are equal but the force on the cube 
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Figure 3. The lateral energy for two gold square-based prisms. 

is only half the size. The result for a cylinder with height L and a circular base with 
radius R = L/y/n follows closely the result for the cube, with a deviation smaller than 
the numerical error ~ 1%. These results show that the shape of the side that faces 
the substrate is not important. One consequence of this result is the additivity of the 
vdW energy for prisms: At all distances the vdW energy for prisms of same height 
and arbitrary shape is approximately proportional to the area of its base. 

3.2.2. Interactions between two objects: The dependence of the vdW interaction on 
the object shapes becomes more complex in the case of the interaction between two 
objects. The interaction tends to modify the position and orientation of the particles 
to minimize the energy. It is possible to consider configurations in which two important 
kinds of interaction arise: lateral forces and rotational forces. 

In figure[3]we show results for lateral forces between two prisms of height L /2 with 
square bases of side L; they are placed above each other at a fixed distance d; they are 
displaced in the horizontal plane a distance I from the complete parallel alignment. 
The solid line shows the results when the prisms are displaced in the direction parallel 
to one of its side faces. The dashed line shows the result when the prisms are displaced 
in the diagonal direction relative its base (see the lower inset of figure [3l) We note 
that both curves reach a minimum in the position of complete alignment. From this 
one can expect an oscillatory behaviour around this position. Even thought one could 
expect the additivity of the interaction in the limit when the particles are close to each 
other, there is no such additivity with the conditions considered here. In the diagonal 
displacement the area of interaction increases differently from that in the other type 
of displacement. However the energy shows an increase similar to that for the other 
displacement. This result suggests that the surface and boundary effects due to the 
interaction between the surface plasmons of the different objects is more important in 
the generation of lateral forces than in the case of an object above a substrate. 

In figure[4]we show results for the rotational energy between two crossed cylinders 
as shown in the inset. The cylinders are rotated relative each other an angle 9 around 
an axis passing through their centers. The angle varies from to n/2. The upper two 
curves are for cylinders with circular cross section with diameter L. The dotted (solid) 
curve is for cylinders of length L (2L). The lower two curves are for cylinders with 
square cross section of side length L. The dashed (dash-dotted) curve is for cylinders 
of length L (2L). In all cases the closest distance, d, is 0.3L. 

We see that in the case of perpendicular cylinders the interaction energy for short 
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Figure 4. The rotation energy for finite gold cylinders the distance d = 0.3L 
apart. 

and long cylinders are not the same, in disagreement with the prediction from the 
PFA model. This fact shows that the interaction between the parts of the cylinder 
that do not face parts of the other one has a strong contribution to the interaction 
energy. In spite of the fact that the interaction between the circular cylinders shows 
a change of behaviour from that between two spheres, ~ c? _1 , when they are fully 
crossed to that between two parallel cylinders, <~ g?~ 3 / 2 , when they are aligned, we 
observe the strongest changes of energy in the case of square cylinders; this is due to 
the fact that the plane surfaces have a stronger interaction than the curved surfaces. 
We can also see that there is no perfect additivity in this configuration. In the case of 
complete alignment (9 = 0) the cylinders of double length do not show double energy. 
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